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Abstract. Elliptic functions considered by Dixon in the nineteenth century 
and related to Fermat's cubic, + = 1, lead to a new set of continued frac- 
tion expansions with sextic numerators and cubic denominators. The functions 
and the fractions are pregnant with interesting combinatorics, including a spe- 
cial Polya urn, a continuous-time branching process of the Yule type, as well as 
permutations satisfying various constraints that involve either parity of levels 
of elements or a repetitive pattern of order three. The combinatorial models 
are related to but different from models of elliptic functions earlier introduced 
by Viennot, Flajolet, Dumont, and Frangon. 



In 1978, Apery announced an amazing discovery: "C(3) = l/"-"^ irrational" . 
This represents a great piece of Eulerian mathematics of which van der Poorten has 
written a particularly vivid account in [58) . At the time of Apery's result, nothing 
was known about the arithmetic nature of the zeta values at odd integers, and not 
unnaturally his theorem triggered interest in a whole range of problems that are 
now recognized to relate to much "deep" mathematics |2H1 Apery's original 
irrationality proof crucially depends on a continued fraction representation of C(3)- 
To wit: 
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where n7(n) = (2n + l)(17n(n + 1) -I- 5). 

(This is not of a form usually considered by number theorists.) What is of special 
interest to us is that the nth stage of the fraction involves the sextic numerator , 
while the corresponding numerator is a cubic polynomial in n. Mention must also 
be made at this stage of a fraction due to Stieltjes (to be later rediscovered and 
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extended by Ramanujan |4| Ch. 12]), namely, 

1 



k=l 
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(2) ^(1) ^ 

a{2) - - 

where a{n) = (2n + l)(n(n + 1) + 2z{z + 1) + 1). 

Unfortunately that one seems to have no useful arithmetic content. 

Explicit continued fraction expansions of special functions are really rare |48[I6(J| . 
Amongst the very few known, most involve numerators and denominators which, 
at depth n in the continued fraction, depend rationally on n in a manner that is at 
most quartic. In this context, the Stieltjes-Ramanujan- Apery fractions stand out. It 
then came as a surprise that certain functions related to Dixon's nineteenth century 
parametrization of Fermat's cubic + Y^ = 1 lead to continued fractions that 
precisely share with Q and (O the cubic-sextic dependency of the their coefficients 
on the depth n, for instance. 



sm(u)e "/"^ du 
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„ 7 -8^ -9^ -10x6 
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where 6„ = 2(3?i + l)((3n + 1)^ + 1). 

There, the function sm is in essence the inverse of a 2^'i-hypergeometric of type 
|; |]; see below for proper definitions. This discovery, accompanied by several 
related continued fractions, was first reported in Conrad's PhD thesis [E! defended 
in 2002. It startled the second author with a long standing interest in continued 
fractions |19[ I20L 1^. when he discovered from reading Conrad's thesis in early 
2005, that certain elliptic functions could precisely lead to a cubic-sextic fraction. 
This paper describes our ensuing exchanges. We propose to show that there is an 
interesting orbit of ideas and results surrounding the Fermat cubic, Dixon's elliptic 
functions, Conrad's fractions, as well as certain urn models of probability theory, 
Flajolct's combinatorial theory of continued fractions, and finally, the elementary 
combinatorics of permutations. 

Plan of the paper. The Dixonian elliptic functions, "sm" and "cm, are intro- 
duced in Section^ Their basic properties are derived from a fundamental nonlinear 
differential system that they satisfy. It this way, one can prove simply that they 
parametrize the Fermat cubic and at the same time admit of representations as 
inverses of hypergeometric functions. Section 13 presents the complete proof of six 
continued fractions of the Jacobi type and three fractions of the Stieltjes type that 
are associated (via a Laplace transform) to Dixonian functions — this is the first 
appearance in print of results from Conrad's PhD thesis jj^l- Next, in Section O 
we prove that the combinatorics of the nonlinear differential system defining sm, cm 
is isomorphic to the stochastic evolution of a special process, which is an urn of the 
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Polya type. As a consequence, this urn, together with its continuous-time analogue 
(a two-particle version of a classical binary branching process) can be solved ana- 
lytically in terms of Dixonian functions. Conversely, this isomorphism constitutes 
a first interpretation of the coeSicients of sm, cm phrased in terms of combina- 
torial objects that are urn histories. Our second combinatorial interpretation is 
in terms of permutations and it involves peaks, valleys, double rises, and double 
falls as well as the parity of levels of nodes in an associated tree representation; 
see Section 0] (Several combinatorial models of elliptic functions due to Dumont, 
Flajolet, Frangon, and Viennot I14L 1151 ITOl 1211 159| are otherwise known to involve 
parity-constrained permutations.) Next, the continued fraction expansions relative 
to Dixonian functions can be read combinatorially through the glasses of a theorem 
of Flajolet and a bijection due to Frangon and Viennot relative to systems 
of weighted lattice paths and continued fractions. This is done in Section which 
presents our third combinatorial model of Dixonian functions: the coefficients of 
sm, cm are shown to enumerate certain types of permutations involving a repetitive 
pattern of order three. Finally, Section El briefly summarizes a few other works 
where Dixonian functions make an appearance. 

Warning. This paper corresponds to an invited lecture at the Viennotfest and, 
as such, its style is often informal. Given the scarcity of the literature relative to 
Dixonian functions, we have attempted to provide pointers to all of the relevant 
works available to us. Thus, our article attempts to kill three birds with one stone, 
namely be a tribute to Viennot, survey the area, and present original results. 

Dedication. This paper is kindly dedicated to Gerard. . .Xavier Viennot on the 
occasion of his sixtieth birthday, celebrated at Lucelle in April 2005. His works in 
lattice path enumeration, bijective combinatorics, and the combinatorics of Jaco- 
bian elliptic functions have inspired us throughout the present work. 

1. The Fermat curves, the circle, and the cubic 
The Fermat curve F™ is the complex algebraic curve defined by the equation 

+ = 1. 

(Fermat- Wiles asserts that this curve has no nontrivial rational point as soon as 
m > 3.) 

Let's start with the innocuous looking F2, that is, the circle. Of interest for this 
discussion is the fact that the circle can be parametrized by trigonometric functions. 
Consider the two functions from C to C defined by the linear differential system, 

(4) s' = c, c' = — s, with initial conditions s(0) = 0, c(0) = 1. 

Then the transcendental functions s, c do parametrize the circle, since 

s{z)^+c{zf = 1, 

as is verified immediately from the differential system, which implies (s^ -l-c^)' = 0. 
At this point one can switch to conventional notations and set 

s{z) = smz, c{z) = cos{z). 
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Figure 1. The two Fermat curves F2 and F3. 



It is of interest to note that these functions are also obtained by inversion from an 
Abelian integral^ on the F2 curve: 

/■sin z _ 

TT^F' ^^^^^^ -sin(z)2 
For combinatoriahsts, it is of special interest to note that the related functions 

sin z 1 

tan z = , sec z = , 

cos z cos z 

enumerate a special class of permutations, the alternating (also known as "up-and- 
down" or "zig-zag") ones. This last fact is a classic result of combinatorial analysis 
discovered by Desire Andre around 1880. 

1.1. The Fermat cubic and its Dixonian parametrization. Next to the circle, 
in order of complexity, comes the Fermat cubic F3. Things should be less elementary 
since the Fermat curve has (topological) genus 1 , but this very fact points to strong 
connections with elliptic functions. 

The starting point is a clever generalization of l@J. Consider now the nonlinear 
differential system 

(5) s' = c2, c' = 

with initial conditions s(0) = 0, c(0) — 1. These functions are analytic about the 
origin, a fact resulting from the standard existence theorem for solutions of ordinary 
differential equations. Then, a one line calculation similar to the trig function case 
shows that 

s{zf + c{zf = 1, 

since 

(s3 + c3)' = 3,s2c2-3c2s2 = 0. 
Consequently, the pair (s(z), c{z)) parametrizes the Fermat curve F3, at least locally 
near the point (0, 1). The basic properties of these functions have been elicited by 



^Given an algebraic curve P{z, y) = 0, an Abelian integral is any integral J R{z, y) dz, where 
ij is a rational function. 




Alfred Cardew Dixon (1865-1936) in a long paper [T^. Dixon established that the 
functions are meromorphic in the whole of the complex plane and doubly periodic 
(that is, elliptic), hence they provide a global parametrization of the Fermat cubic. 

From now on, we'll give the s, c functions the name introduced by Dixon (see 
Figure El for a rendering) and set 

sm(z) = s(z), cm(z) = c(z). 

Their Taylor expansions at (not currently found in Sloane's Encyclopedia of In- 
teger Sequences (EIS) |54j ) start as follows: 

{^4 ^7 ^10 ^13 
sm(z) = z - 4— + 160— - 20800— + 6476800— 
cm(z) = l_2-^+40^- 3680^+8880000|^---- . 

In summary, we shall operate with the new notations sm, cm and with the defining 
system: 

(7) sm'(z) = cm^(z), cm'(z) = - sm^(z); sm(0) = 0, cm(0) = 1. 

1.2. A hypergeometric connection. At this point, it is worth noting that one 
can easily make the s = sm and c = cm functions somehow "explicit" Start from 
the defining system (System Q abbreviated here as (E)) and apply differentiation 
id): 

(8) s' = ^ ^ s" = 2cc' ^ s" = -2cs' ^ s" = -2s2V7. 



Then "cleverly" multiply by VV to get, via integration (J), 

(9) s"V7=-2sV ^ + 

3 3 

for some integration constant K which must be equal to | , given the initial condi- 
tions. This proves in two lines(!) that sm is the inverse of an integral, 

(10) X = 
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this integral being at the same time an incomplete Beta integral and an Abclian 
integral over the Fermat curve (/ By the same devices, it is seen that the 
function cm(z) also satisfies 

and is thus itself the inverse of an Abelian integral. 

Via expansion and term-wise integration, the latest finding ()10|l can then be 
rewritten in terms of the classical hypergeometric function, 

^. ^ , a-l3z a(a + 1) • /3(/3 + 1) z2 

7 1! 7(7 + 1) 2! 

with special rational values of the parameters. Letting Inv(/) denote the inverse 
of / with respect to composition (i.e., Inv(/) — giifog^g°f~ Id) we state: 

Proposition 1 (Dixon 13 ). The function sm is the inverse of an Abelian integral 
over F3 and equivalently the inverse of a 2F1: 

sm(z)=Inv/ =Invz-2F, z 



The function cm is then defined near by cm(z) = s^l — sm^(z), or alternatively 
by (fTHl. 

The analogy with the sine function is striking. Of course. Proposition is not 
new and all this is related to extremely classical material. Dixon jEl discusses the 
implicit integral representations (|10|l . Hll|) and writes concerning the prehistory of 
his sm, cm: 

The only direct references that I have come across elsewhere are certain 
passages in the lectures of Professor Cayley and Mr. Forsyth where 
the integral J — — '''3^2/3 '^^^ used to illustrate Abel's Theorem, in the 
treatises of Legendre, and Briot and Bouquet, and again in Professor 
Cayley's lectures and elsewhere where it is shewn how to turn the integral 
into elliptic form, and lastly at the end of Bobek's Einleitung in die 
Theorie der elliptischen Functionen where expressions are found for the 
coordinates of any point on the above cubic. 

It is fascinating to be able to develop a fair amount of the theory from the differential 
equation ITJ, using only first principles of analysis. (See Dixon's article as well as 
Sectional of the present paper for more.) 

Note 1. Lundberg's hypergoniometric functions. The question of higher degree general- 
izations of the differential system satisfied by sm, cm is a natural one. Indeed, the system 
s' = c''~^ ,c' = —s^~^ parametrizes locally near (0, 1) the Fermat curve X'' + Y'' — 1. (This 
parametrization ceases to be a global one, however, since Fp has genus g — (p— l)(p — 2)/2, 
so that <; > 3 as soon as p > 4; see for instance |39| . Thus, only the case p = 3 leads directly 
to elliptic functions.) The corresponding functions are still locally inverses of Abelian in- 
tegrals over the Fermat curve, which is verified by calculations similar to Equations (|5J 
and @. In fact, a rather unknown high school teacher from Sweden, Erik Lundberg, 
developed in 1879 an elementary theory (see |44p of what he called "hypergoniometric 
functions", including a sinualis and a cosinualis that are indexed by rational numbers. 
Interest in these questions was recently rekindled by an insightful article of Lindqvist and 
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Peetre published in the American Mathematical Monthly |42) . The authors discuss vari- 
ous connections to eUiptic functions, in particular, the reduction of Fermat's cubic to its 
Weierstrafi normal form by elementary manipulations. See also a problem 4r in the same 
issue of the Monthly posed by these authors together with a remark due to Jon Borwein 
to the effect that 

sm ("^ ^ ~ cm(a::). 
There 713 is a fundamental constant of Dixonian functions: 

Many more properties of this type are to be found in Dixon's paper. In particular, one 
has, 

sm(7r3 + u) = sm(it), cm(7r3 + u) — cm{u), 
that is, TTs is a real period. The complete lattice of periods of sm, cm is 

which is consistent with the rotational invariance: sm{ijju) — lu sm(ii) and cm(a;u) = 
cm(M). It thus corresponds to the hexagonal lattice displayed on Figure |S| and to case C 
of the urn models evoked in Subsection 13.51 □ 



2. Some startling fractions 

In 1907, L. J. Rogers devised two methods to obtain continued fraction 
expansions of Laplace transforms of the Jacobian elliptic functions sn, cn (see for 
instance [HJj for the basic theory of these functions). In his first method, he re- 
sorted to integration by parts. His second method involves a general "addition 
theorem" 03 ; it was to some extent a rediscovery of a technique introduced 
by T. J. Stieltjes in j56| that relies on diagonalization of certain infinite quadratic 
forms. Stieltjes and Rogers found altogether three families of continued fractions 
relative to sn, cn. (A fourth family, implicit in the work of Rogers, was discovered 
almost a century later by Ismail and Masson 33 .) Such expansions turn out to 
be useful: S. Milne in '46' obtained additional expansions implying explicit Hankel 
determinant evaluations, which enabled him to prove some deep results about sums 
of squares and sums of triangular numbers. Milne's results in particular include 
exact explicit infinite families of identities expressing the number of ways to write 
an integer as the sum of 4A^^ or 4:N{N + 1) squares of integers, where N is an 
arbitrary positive integer. 

In his PhD thesis defended in 2002, Conrad !121 Ch. 3] applied the integration by 
parts method to develop completely new continued fractions arising from Dixon's 
elliptic functions. These fractions fall into six families which are naturally grouped 
as two sets of three. The underlying symmetries suggest that these classifications 
are fundamental and complete. 

In what follows, we make use of the Laplace transform classically defined by 

C{f,s)^ / /(u)e-^"du. 
Jo 

For our purposes, it turns out to be convenient to set s — x~^. In that case, one 
has 
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which means that the Laplace transform formaUy maps exponential generating 
functions (EGFs) to ordinary generating functions (OGFs): 



(13) -'^^|E«^^'--M-E 



There are a number of conditions ensuring the analytic or asymptotic validity of 
this last equation. In what follows, we shall make use of the integral notation for 
the Laplace transformation, but only use it as a convenient way to represent the 
formal transformation from EGFs to OGFs in (|13|l . 

The continued fractions that we derive are of two types |48[ I60| . The first type, 
called a J-fraction (J stands for Jacobi), associates to a formal power series f{z) a 
fraction whose denominators are linear functions of the variable z and numerators 
are quadratic monomials: 

/W = ^ 



biz' 

1- CqZ 



U2i 



1 — ClZ ■ 

(Such fractions are the ones naturally associated to orthogonal polynomials '11'.) 
The second type^, called in this paper an S'-fraction (5 stands for Stieltjes) has 
denominators that are the constant 1 and numerators that are monomials of the 
first degree: 



1 



1 



diz 



I 'hi. 



An S'-fraction can always be contracted into a J-fraction, with the corresponding 
formulae being explicit, but not conversely. Accordingly, from the point of view of 
the theory of special functions, an explicit S'-fraction expansion should be regarded 
as a stronger form than its J-fraction counterparts. 

2.1. J-fractions for the Dixon functions. We introduce three classes of formal 
Laplace transforms: 



S„ := / sm(u)"e""/^dw 

(14) Cn := / sm(u)" cm(u) e^"/^dw 

Jo 

/•oo 

Dn ■■= / sm(u)"cm(u)2e-"/^du 







(Since sm'^ -I- cm'^ = 1, we can reduce any polynomial in sm and cm to one which is 
at most of degree 2 in cm. We let powers of sm grow since smu vanishes at u = 0.) 



^An alternative name for S'-fraction is "regular C-fraction" . What we called J-fraction is also 
known as an "associated continued fraction". 
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After integrating by parts and reducing to canonical form using the Fermat cubic 
parametrization, we obtain the following recurrences: 

5*0 = X, Sn = nxDn^i (n > 0) 

(15) Co = x-xS2, Cn ^ rixSn-i - {n + l)xSn+2 (n > 0) 

Do = X — 2xG2, Dn = nxCn-i — {n + 2)xCn+2 {n > 0) 

These recurrences will provide six J- fractions and three iS-fr actions. 

For Laplace transforms of powers of sm(u), the recurrences (I15II yield the fun- 
damental relation, 

Sn n{n - l){n - 2)x^ 



(16) 



Sn-3 l + 2n{n^ + l)x^ -n{n+l){n + 2)x^^^ 
subject to the initial conditions 

x^ 2x^ _ 6x^ 

'-'1 ^ Z ; — n „ □ c. ! >^2 — , — — "33 — 



1 + 4x3 - 6a;3|^' 1 + 20a;3 - 24a;3 1^ ' 1 + 60a;3 - 60x3 1| ' 

Repeated use of the relation (|16|l starting from the initial conditions then leads to 
three continued fraction expansions relative to each of Si, 82, S^. 

An entirely similar process applies to C„. We find 

Cn _ n(n- l)(n- 2)a;3 

Cn-3 ~ 1 + ((n - l)n2 + {n + l)2(n + 2)) x^ - {n + l)(n + 2)(n + 3)^3 ' 

subject to 

^ Q ^ 

°~ 1 + 2x3-6x3^' ^ ~ 1 + 12x3 - 24x3^' ^ ~ i + ^q^s ^ qq^^C^' 
Oo Oi 02 

For Dn, we could handle things in a similar manner, but it turns out that there 
is a simpler procedure available to us. Differentiating powers of sm(u) with respect 
to It, we have 

sm(u)"+^ ^ (n + 1) sm(u)" cm^(u). 

du 

Given the effect of Laplace transforms on derivatives, it follows that 

/:(sm(u)" cm(u)2,x-i) = ^-—C{sm''+\u),x^^). 

^ ^ ' ^ ' ' (n + l)x 

We thus obtain no substantially new J-fractions from D„ by this process. 

Theorem 1 (Conrad |12|'). The formal Laplace transform of the functions sm" 
(n = 1,2,3) and sm" • cm (n = 0,1,2) have explicit J -fraction expansions with 
cubic denominators and sextic numerators as given in Figure\^ 

Consider the list of positive integers in ascending order with each integer listed 
twice: 

1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, . . . 

There are essentially six ways to break this up into ascending 6-tuples, allowing for 
the possibility of missing leading entries in the first 6-tuple. The a„ terms in the J- 
fraction for the Laplace transform of sm(u), sm^{u), sm3(u) are seen to correspond 
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POO 

1 SVO.U e""'''^ du = 
Jo 








aix^ 

a2X 

1 + bix^ 

1 + 62X3 - ■ ■ . 


where 


On 


= (3n - 2)(3n - l)^{3nf{3n 


+ 1), 6„ = 2(3n + l)((3n + 1)2 + 1). 






poo 

/ sm2(w)e-"/^dM = 


2x^ 






aix^ 

1 + 20x^ 

aix^ 

1 + 6ia;3 

1 + 62x3 - ■ • ■ 


where 


ttn 


= (3n - l)(3n)^(3n + lf{3n 


+ 2), 6n = 2(3n + 2)((3n + 2)2 + l). 






POO 

/ sm^(w)e-"''^dw = 








OlX® 

1 1 1 

1 -f- yjMX -\ ^ 

1lX 

1 + 6ix3 

1 + 62X3 - ■ • . 


where 




= (3n)(3n + lf{3n + 2f{3n 


+ 3), 6„ = 2(3n + 3)((3n + 3)2 + l). 






POO 

/ cm(w) e""/^ du = 
Jo 


X 






aix® 

1 -f- ZX -f- 

02X 

1 + 6ix3 

1 + 62X3 - ■ ■ . 


where 


a„ 


= (3n-2)2(3n-l)^(3n)^ 


6„ = (3n - l)(3n)^ + (3n + l)^(3n + 2). 






roo 

/ sm(w) cm(u) e""''^ du = 
Jo 


x^ 






OlX® 

1 1 1 '^t'^ 

a2X 

1 + 6ix3 

1 + 62x3- 


where 


an 


= (3n- lf{3nf{3n+lf, 


6„ = (3n)(3n + 1)^ + (3n + 2f(3n + 3). 






f oc 

/ sm^(M) cm(w) e~"^^ du = 
Jo 


2x^ 






OlX® 

1 + 40x3 + g 

02X 

1 + 6ix3 

1 + 62X3 - ■ • . 


where 




= {3nf{3n + lf{3n + 2f, 


6„ = (3n + l)(3n + 2f + (3n + 3)^(371 + 4). 



Figure 3. The six basic J- fractions relative to Dixonian functions. 
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to the partitions: 

sm{u) : (1, 2, 2, 3, 3, 4), (4, 5, 5, 6, 6, 7), (7, 8, 8, 9, 9, 10), . . . , 
sni2(^.) : (2, 3, 3, 4, 4, 5), (5, 6, 6, 7, 7, 8), (8, 9, 9, 10, 10, 11), . . . , 
sni3(^.) : (3, 4, 4, 5, 5, 6), (6, 7, 7, 8, 8, 9), (9, 10, 10, 11, 11, 12), . . . , 

These are the three possible "odd" partitions, odd in the sense that the first integer 
in each 6-tuple appears exactly once. The three remaining continued fractions, those 
of cm, cm • sm, cm • sm^ , are associated in the same way with the three possible even 
partitions. 



2.2. S'-fractions for Dixon functions. Starting with the original recurrences of 
the previous section(Equation ifTCjl . we can use the relation between Sn and Dn^i to 
eliminate either the letter S or the letter D. These recurrences reduce to S'- fraction 
recurrences, which we tabulate here: 

Sn n{n ~ l)x^ Cn nx 

Cn-2 ^ l + n(n+l)a;2%±i' S'„-i ^ 1 + (n + l)a;%^ 



Dn-2 i + („ + i)(„ + 2)a;2%±i' C„_i l + (n + 2)a;^' 



+2 



We need initial conditions to generate S'-fractions, and the six candidate starting 
points give just three 5-fraction initial conditions: 

Si = x^-2x^C2 



01 



Co = X — xS; 



Ci = x^- 2xS3 



l + 2x^ 



The three that fail to give good initial conditions for an S"- fraction are as follows: 

2x^ 



52 = 2x^ - 2x^S2 - 6x2(73 

53 = 6a;4 - 12x^S3 - Ux^d 
C2 = 2x^ - Ax^C2 - ix'^Si 



l + 2a;3 + 6a;2^ 
Qx^ 

1 + 12.i;3 + 12x2^ 
2x^ 

1 + 4x3 + 3x2^' 



On iterating the first three relations, we obtain three S-fraction expansions: 

Theorem 2 (Conrad jl2|). The formal Laplace transform of the functions sm, cm 
and sm • cm have explicit S-fraction expansions with cubic numerators, 

smue ^ du = 



aix^ 
1 + — 



1 + ^ 
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where for r > 1; aar-i = (3r - 2)(3r - l)^, a2r = {3rf{3r + 1); 

cmue "'^du^ ^ 

aix 
1 + 



1 + ^ 



1+ ■■• 

where for r> I: as^-i = (3r - 2)2(3r - 1), = (3r - l)(3r)2; 



oo 

sm 



(m) cm(u) e "/"^ 



1 + — 



1 + ^ 



1+ 

where for r>\: asr-i = (3r - l)2(3r), = (3r)(3r + 1)2. 



Consider now the set of positive integers written in increasing order with each 
integer written twice: 

1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, . . . 

There are three ways to divide this hst into ascending triples if we permit ourselves 
to discard leading ones, corresponding to the numerators coefficients a„ of the S- 
fractions associated to cm, sm, and sm • cm: 

cm(u) : (1,1,2), (2,3,3), (4,4,5), (5,6,6), ... , 

sm{u) : (1,2,2), (3,3,4), (4,5,5), (6,6,7), ... , 

sm{u) ■ cm(u) : (2, 2, 3), (3, 4, 4), (5, 5, 6), (6, 7, 7), ... . 

This correspondence suggests that these three continued fractions form a complete 
set. 



Note 2. The cubic + - 3aXY = 1. Conrad in his dissertation follows Dixon 
and examines the larger class of functions corresponding to the cubic 

+ - 3aXY = 1, 

for arbitrary a. These do give rise to continued fractions, but ones that are non-standard: 
they are not of the or J type as they involve denominators that are linear (in x) and 
numerators that are cubic. 

□ 



An alternative derivation of the J-fraction expansions, based on a direct use of 
the differential system, will be given in Subsection 15. II when we discuss a method 
of Andre (Note|Hl). 
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3. First combinatorial model: Balls games 

Dixonian elliptic functions, as we'll soon see, serve to describe the evolution 
of a simple urn model with balls of two colours. In this perspective, the Taylor 
coefficients of sm, cm count certain combinatorial objects that are urn histories, 
or equivalently, weighted knight's walks in the discrete plane — this provides our 
first combinatorial interpretation. Going from the discrete to the continuous, we 
furthermore find that Dixonian functions quantify the extreme behaviour of a classic 
continuous-time branching process, the Yule process. Finally, we show that the 
composition of the system at any instant, whether in the discrete or the continuous, 
case, can be fully worked out and is once more expressible in terms of Dixonian 
functions. 

3.1. Urn models. Balls games have been of interest to probabilists since the dawn 
of time. For instance in his Theorie analytique des probabilites (first published 
in 1812), Laplace writes "C/ne urne A renfermant un tres grand nomhre n de boules 
blanches et noires; a chaque tirage, on en extrait une que Von remplace par une boule 
noire; on demande la prohahilite qu'apres r tirages, le nomhre des boules blanches 
sera x." Such games were later systematically studied by Polya. What is directly 
relevant to us here is the following version known in the standard probability liter- 
ature as the Polya urn model (also Polya-Eggenberger) : 

Polya urn model. An urn is given that contains black and white balls. 
At each epoch, a ball in the urn is chosen at random (but not removed). If 
it is black, then a black and [3 white halls are placed into the urn; else it is 
white and 7 black and S white halls are placed into the urn. 

The model is fully described by the "placement matrix" , 

The most frequently encountered models are balanced, meaning that a + f3 = j + S, 
and negative entries in a matrix M are interpreted as subtraction (rather than ad- 
dition) of balls. For instance, Laplace's original problem corresponds to -^^l)' 
which is none other than the coupon collector's problem in modern terminology. 
What is sought in various areas of science is some characterization, exact or asymp- 
totic, of the composition of the urn at epoch n, given fixed initial conditions. The 
elementary introduction by Johnson and Kotz |36| mentions applications to sam- 
pling statistics, learning processes, decision theory, and genetics. Recently, Polya 
urn models have been found to be of interest in the analysis of several algorithms 
and data structures of computer science; see especially Mahmoud's survey 

The main character of this section is the special urn defined by the matrix 

M12 = 

This can be visualized as a game with balls of either black ('x') or white ('y') 
colour. If a ball is chosen, it is removed [the —1 entry in the matrix] from the urn 
and replaced by two balls of the opposite colour according to the rule 

(17) X — > yy, y — > xx. 
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A history of length n (see Frangon's work 26, for this terminology and similar 
ideas) is, loosely speaking, any description of a legal sequence of n moves of the 
Polya urn. Let conventionally the urn be initialized at time with one black ball (x). 
A history (of length n) is obtained by starting with the one- letter word x at time 
and successively applying {n times) the rules of (I17|l . For instance, 

X — > yy — > yxx — > yyyx — » xxyyx — > xyyyyx, 

is the complete description of a history of length 5. (The replaced letters have been 
underlined for readability.) We let Hn^k be the number of histories that start with 
an X and, after n actions, result in a word having k occurrences of y (hence n+l — k 
occurrences of x) . Clearly the total number of histories of length n satisfies 

k 

since the number of choices is 1, 2, 3, 
all histories of length < 3: 

n — : x 
n — 1 : x 
n — 2 : X 

X 

n = 3 : X 

X 
X 
X 
X 
X 

The sequence (ffn.o) starts as 1,0,0,2 for n — 0,1,2,3 and it is of interest to 
characterize these combinatorial numbers. 



. . . at times 1,2,3.... Here is a small table of 



yy 






yy 


— > xxy 




yy 


— > yxx 




yy 


— > xxy — > 


yyxy 


yy 


— > xxy — > 


xyyy 


yy 


— > xxy — > 


xxxx 


yy 


— > yxx — > 


xxxx 


yy 


— > yxx — > 


yyyx 


yy 


— > yxx — > 


yxyy 



3.2. Urns and Dixonian functions. Let us come back to Dixonian functions. 
Consider for notational convenience the (autonomous, nonlinear) ordinary differen- 
tial system 

(18) S: —=y^^ with x{Q) = xo, y{0) = yo, 

which is the signless version of ((^l. In this subsection, we only need the spe- 
cialization a;(0) = 0, y{0) = 1, but we will make use of the general case H18|l in 
Subsection 13.41 below. 

The pair {x{t),y{t)) with initial conditions .t(0) = 0, j/(0) = 1 parametrizes the 
"Fermat hyperbola" , 

- = 1, 

which is plainly obtained from Fermat 's "circle" F3 by a vertical symmetry. It 
is natural to denote this pair of functions by (smh(t), cmh(t)), these functions 
smh, cmh being trivial variants of sm, cm, where the alternation of signs has been 
suppressed: 



(19) 



smh(z) — — sm(— 2:), cmh(z) = cm(— z). 
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A first combinatorial approach to Dixonian functions can be developed in a 
straightforward manner by simply looking at the basic algebraic relations induced 
by the system S of (|18|) . To this purpose, we define a linear transformation 5 
acting on the vector space C[x,y] of polynomials in two formal variables x,y that 
is specified by the rules, 

(20) S[x]^y'^, S[y]^x^, S[u ■ v] ^ S[u] ■ v + u ■ S[v], 

u,v being arbitrary elements of C[a;, y]. A purely mechanical way to visualize the 
operation of 5 comes from regarding 5 as a. rewriting system 

(21) X yy, y xx, 

according to the first two rules of (|20|l . The third rule means that (S is a derivation, 
and it can be read algorithmically as follows: when b is to be applied to a mono- 
mial w of total degree d in x, y, first arrange d copies of where in each copy one 
instance of a variable is marked (underlined), then apply the rewrite rule (|21l) once 
in each case to the marked variable, and finally collect the results. For instance, 

xyy i-> xyy, xyy, xyy {yy)yy, x{xx)y, xy{xx), so that 6[xy^] ^ y'^ + 2x^y. 

Two facts should now be clear from the description of 6 and the definition of 
histories. 

(i) Combinatorially, the nth iterate S"[x°'y''] describes the collection of all the 
possible histories at time n of the Polya urn with matrix A^i2, when the 
initial configuration of the urn has a balls of the first type (x) and b balls 
of the second type (y). This is exactly the meaning of the replacement 
rule H21|) . In particular, the coefficient^, 

H^^^f = [x'^y'] 5" [x'^y''] , k + l^n + a + b 

is the number of histories of a Polya urn that lead from the initial state 
x°-y'' to the final state x^y^. 
(ii) Algebraically, the operator S does nothing but describe the "logical conse- 
quences" of the differential system S. In effect, the first two rules of H21() 
mimic the effect of a derivation applied to terms containing x = x(t) and 
y — y{t), "knowing" that x' — y^, y' = x^ . Accordingly, the quantity 
5'^[x°'y^] represents an nth derivative, 

6^[x%f]^^x{try{tf expressed in x{t),y{t), 

where x{t),y{t) solve the differential system x' = y^, y' — x^ . 

In summary, we have a principle^ whose informal statement is as follows. 

Equivalence principle. The algebra of a nonlinear autonomous system 
that is monomial and homogeneous, like S in p8|l . is isomorphic to the 
combinatorics of an associated Polya urn. 

x^j/", then the notation is used to represent coefficient extraction: 

■^In this paper, we have chosen to develop a calculus geared to Dixonian functions. Our 
approach is in fact appreciably more general and it can be applied to any balanced urn model, as 
considered in 1221 . 
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This principle can be used in either direction. For us it makes it possible to analyse 
the Polya urn Aii2 in terms of the functions sm, cm that have already been made 
explicit in Proposition ^ In so doing, we rederive and expand upon an analysis 
given in a recent paper on analytic and probabilistic aspects of urn models 
Here is one of the many consequences of this equivalence. 

Theorem 3 (First combinatorial interpretation) . The exponential generating func- 
tion of histories of the urn with matrix M.12 that start with one hall and terminate 
with balls that are all of the other colour is 

E^»,o^--h(z).^.-sm(-z). 
^-^ nl cm(z) 

The exponential generating function of histories that start with one ball and termi- 
nate with halls that are all of the original colour is 

Hn,n+i—r = cmh(z) = — ■ = cm(-z). 

^ — ' n! cm z 

n>0 ^ ' 

Proof. This is nothing but Taylor's formula. Indeed, for the first equation, we have 
by the combinatorial interpretation of the differential system: 

But from the algebraic interpretation, 

jn 

'5"NL=o..=i= ;^snihW 

The result then follows. □ 

Theorem |31 thus provides a first combinatorial model of Dixonian functions in 
terms of urn histories. 

Note 3. On Dumont. Though ideas come from different sources, there is a striking 
parallel between what we have just presented and some of Dumont's researches in the 
1980's and 1990's. Dumont gives an elegant presentation of Chen grammars in jl8). 
There, he considers chains of general substitution rules on words; such chains are partial 
differential operators in disguise. For instance, our 5 operator is nothing but 

X 2 d o d 
= x +y TT- 
ay ox 

Dumont has shown in J18] that operators of this type can be used to approach a variety of 
questions like rises in permutations. Bell polynomials, increasing trees, parking functions, 
and Polya grammars. The interest of such investigations is also reinforced by consideration 
of the trivariate operator 

d d d 

yzTT- +ZX— +xy — , 
ox ay oz 

itself related to the differential system 

X = yz, y = zx^ z — xy, 

and to elliptic functions of the Jacobian type, which had been formerly researched by 
Dumont and Schett; see |14l 1151 IS^ and Subsection 13.51 below for further comments on 
the operator point of view. □ 
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multiplicity p 



multiplicity q 



Figure 4. The Bousquet-Melou-Petkovsek moves (left) and their 
version with multiplicities (right) for the M12 urn model. 

Note 4. Knight's walks of Bousquet-Melou and Petkovsek. Theorem|3is somehow related 
to other interesting combinatorial objects (Figure [IJ. Say we only look at the possible 
evolutions of the urn, disregarding which particular ball is chosen. An urn with p black 
balls and q white balls may be represented in the cartesian plane by the point P with 
coordinates (p, g). A sequence of moves then defines a polygonal line, PqPi ■ ■ ■ Pn with 
Po ~ (1,0), each move PjPj+i being of type either /? — (—1, +2) [i.e., two steps East, one 
step South] when a black ball is picked up or a; = (2, —1) [i.e., two steps North, one step 
West] otherwise. This defines a random walk on the first quadrant that makes use of two 
types of knight moves on the chessboard. The enumeration of these walks is a nontrivial 
combinatorial problem that has been solved recently by Bousquet-Melou and Petkovsek 
in |7|. They show for instance that the ordinary generating function of walks that start 
at (1,0) and end on the horizontal axis is 



There, ^^'^ = ^o- • -o^ is the ith iterate of ^. These walks have a merit 0: they provide an 
extremely simple example of a linear recurrence with constant coefficients whose generating 
function is highly transcendental (in fact, not even holonomic). 

In order to obtain a complete history from a knight's walk, one has to add some 
supplementary information, namely, which ball is chosen at each stage. This corresponds 
to introducing multiplicative weights. The rule is then as follows. For a walk with a /3-step 
that starts at point (p, q) (a black ball is chosen) the weight is p; dually, for an tj-step the 
weight is q. □ 

In summary, at this stage, we have available three variants of the interpretation 
of the Taylor coefficients of the Dixonian functions sm, cm provided by Theorem |31 
(i) the enumeration of urn histories relative to the A4i2 urn; (ii) the iterates of 
the special operator 6 = x^dy + y^dx in the style of Dumont; {Hi) multiplicatively 
weighted knight walks of the type introduced by Bousquet-Melou and Petkovsek. 

Note 5. A probabilistic consequence: extreme large deviations. From a probabilistic 
standpoint, the Polya urn model is a discrete time Markov chain with a denumerable set 
of states embedded in Z x Z. The number of black balls at time n then becomes a random 
variable, Xn. Theorem |21 quantifies the probability (P) of extreme large deviations of X„ 



G(x) = ^(-l)"(C«(x)C<'+^>(:c)) 



i>0 

where ^ is a branch of the (genus 0) cubic x^-x^ -(■^ = 0: 
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as 

P(X„ = 0) = ^ = [z"] smh(2) = [z"] - sm(-z). 
m 

Then by an easy analysis ol singularities, one finds that this quantity decreases exponen- 
tially fast: 

in accordance with p2|l . For instance, we have to lOD: 

1/3 

= 1.76663 87502 ■ ■ ■ ; y = 1.76663 87490 ■ • • . 

One can refer to the calculation of 22 , but is is just as easy to note that Proposition Q 
provides directly the dominant positive singularity p of smh(z) as a special value of the 
fundamental Abelian integral. □ 



[z""] sm(z) 
[z^^] sm{z) 



3.3. A continuous-time branching process. The Dixonian functions also make 
it possible to answer questions concerning chain reactions in a certain form of 
particle physics. You have two types of particles, say, foatons and viennons. Any 
particle lives an amount of time T that is an exponentially distributed random 
variable (i.e., P(r > t) — e^*), this independently of the other particles; then 
it disintegrates into two particles of the other type. Thus a foaton gives rise to 
two viennons and a viennon gives rise to two foatons (see Figure ISJ. What is the 
composition of the system at some time t > 0, assuming one starts with one foaton? 

The resemblance with the Polya urn should be clear (with a foaton being a black x 
and a viennon a white y). First, considering only the total population of particles, 
one gets what is perhaps the simplest continuous-time branching process, classically 
known as the Yule process. (For useful information regarding this process, see for 
instance the recent articles by Chauvin, Rouault, and collaborators (HIEII-) Let 
Sk{t) be the probability that the total population at time t is of size k. Introduce 




Figure 5. A truncated view of a tree of particles provided by the 
Yule process (the ordinate of each particle represents the time at 
which it splits into two particles of the other type) . 
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the bivariate generating function, 

oo 

Examination of what happens between times and dt leads to the recurrence 
(22) Sk{t + dt) = (1 - dt)Sk{t) +dt S^{t)SJ{t), 

i+j=k 

which is none other than the usual "backwards equation" of Markov processes. 
Then, one has 

Slit) + Skit) ^ V s.,it)s,it). 



E 

i+j=k 



This induces a nonlinear equation satisfied by S, namely, 

E'it; w) + S(t; w) = S(t; ^^;)^ S(0, w) = w, 

where derivatives are implicitly taken with respect to the time parameter t. The 
solution of this ordinary differential equation is easily found by separation of vari- 
ables, 

X we"* 
1 — wil — e 

which yields 

(23) Skit) ^ e-' il ~ e-')"-', fc > 1. 

In summary, the size of the population at time t obeys a geometric law of parameter 
(1 — e~*), with expectation e*. (The previous calculations are of course extremely 
classical: they are disposed of in just six lines in Athreya and Ney's treatise 
p. 109].) 

The result of H23|l shows that any calculation under the discrete urn model can 
be automatically transferred to the continuous branching process. Precisely, let Ti. 
be the set of all histories of the Polya urn with matrix A4i2- Let /C C 7i be a subset 
of Ti. and let Kiz) be the exponential generating function of /C. Then, by virtue 
of (|23() . the probability that, at time t, the Yule process has evolved according to a 
history that lies in K. is given by 

(24) e-*X(l-e"*). 
This has an immediate consequence. 

Proposition 2. Consider the Yule process with two types of particles. The proba- 
bilities that particles are all of the second type at time t are 

Xit) = e"* smh(l - e"*), y(t) ^ e"* cmh(l - e"'), 

depending on whether the system at time is initialized with one particle of the 
first type (X) or of the second type (Y ). 

Like its discrete-time counterpart, this proposition can be used to quantify ex- 
treme deviations: the probability that all particles be of the second type at time t 
is asymptotic to 

(25) e"*smh(l), smh(l) = 1.20541 51514. 
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At time t, the Yule system has expected size e*. On the other hand, for the discrete 
process, the probabihty at a large discrete time that balls or particles be all of 
one colour is exponentially small. These two facts might lead to expect that, in 
continuous time, the probability of all particles to be of the same colour is doubly 
(and not singly, as in Equation pSfl l exponentially small. The apparent paradox 
is resolved by observing that the extreme large deviation regime is driven by the 
very few cases where the system consists of 0(1) particles only, an event whose 
probability is exponentially small (but not a double exponential). 

Note 6. An alternative direct derivation. The memoryless nature of the process implies, 
like in (I22L the differential system: 

X{t) + X'{t) = Y{tf, Y{t) + Y'{t) = Xitf, 

with initial conditions X{0) = 0, Y{0) = 1. This system closely resembles some of our 
earlier equations. By simple algebra [first multiply by e*, then apply the change of variables 
t l—^ e~*], the solution in terms of Dixonian functions results. □ 



3.4. Dynamics of the Polya and Yule processes. The previous sections have 
quantified the extreme behaviour of the processes — what is the probability that 
balls be all of one colour? In fact, a slight modification of previous arguments give 
complete access to the composition of the system at any instant. 

Consider once more the differential system 

dx dij 

(26) -dt^y^' ^'^^ a;(0) = a;o, y{0) = yo 

(this repeats Equation l|18|) '). The initial conditions are now treated as free param- 
eters or formal variables. This system can be solved exactly by means of previ- 
ously exposed techniques, and, by virtue of the equivalence principle, its solution 
describes combinatorially the composition of the urn, not just the extremal config- 
urations (that correspond to xq = 0, j/o = !)• 

First, the solution to the system. Following the chain of JHJ and ©, one finds 

x'(t)3/2 ^ ^(^)3 ^ ^3^ ^3 ^3 _ ^3^ 

where the initial conditions have been taken into account. This last equation can 
be put under the form 



which integrates to give 
(27) At = 



X3 + A3)2/3 

^/'^ dw 



1, 



The function smh(i) is also defined by an inversion of an integral, 

^ dw 

(1 + ^3)2/3' 

as is easily verified by the technique of Section H"^ Then Equation 1)27(1 provides 
(28) x{t)^Asmh(At + l(^ 



t — /(smh(t)), where I{y) := / 

^0 



A 



which constitutes our main equation. 
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Figure 6. A plot of the coefficients of ;jri5"[a;] for n = 50 and 
n = 150 illustrating the asymptotically gaussian character of coef- 
ficients. 

Second, the interpretation of S in Section 13.21 implies that (5" [x] expresses the 
composition of the urn at the nth stage of its operation when it has been initialized 
with one ball of the first type (x) . Thus the quantity 

j-n 

(29) F(i,xo)-5]^(5"N)_„,^^, 

n>0 

is none other than the x{t) solution to the system ()26|l initialized with a;(0) = xq 
and y{0) — 1. On the other hand, 6^[x] represents the nth derivative of x = x{t) 
expressed as a function of x and y. 

There results from Equations H28|l . (|29|l and the accompanying remarks an ex- 
pression for the bivariate generating function of histories with t marking length 
and Xq marking the number of balls of the first kind. Switching to more orthodox 
notations {xq '—^ x), we state: 

Proposition 3. The composition of the M.12 urn at all discrete instants is described 
by 

= (^x + y^^ + 2x^y^ + {ixy^ + 2x*)^ + ---^ , 

where Hn^k is the number of histories of length n of the urn initialized with x and 
terminating with x'=y"+i"''. 

This statement is the incarnation of Theorem 1 of [221 in the case of the Aii2 urn 
model. From it, one can for instance deduce by an application of the Quasi-powers 
Theorem of analytic combinatorics |22l 1241 m] : The distribution of the number of 
balls of the first kind, equivalently the sequence of coefficients of the homogeneous 
polynomial S"'[x], is asymptotically normal. See Figure El 



22 



ERIC VAN FOSSEN CONRAD AND PHILIPPE FLAJOLET 





(M) 


(S) 




Type 


A 


if-l) 


/ 

X 


— X 


V, y' = x^-' 




also Dixonian 


B 




x' 




y = x?y~^ 


y dx y dy 




C 




x' 


= y' 


y' = X? 


2/ ^ + * ^ 


Dixonian 


D 




x' 


= y' 


y = 




lemniscatic \A'A\ 


E 
F 


(rj) 


x' 

X 


-y' 
= y' 


y = x^y-'^ 

y' = x^y"^ 


^ dx dy 
y dx y Qy 





Figure 7. The six elliptic cases of a 2 x 2 Polya urn: matrix 
[M), system (E), operator ((5), and type, following [22]. 



The transfer from discrete to continuous time afforded by H24|l now permits us 
to deduce the composition of the Yule process from the Polya urn via the transfor- 
mation K{z) I— > e^*^K{\ 



Proposition 4. In the Yule process with two types of particles, the probability 
generating function of the number of particles of the first type (x ) at time t is 



e 



(l + s3)2/3 



By means of the continuity theorem for characteristic functions, it can be verified 
that the distribution of the number of particles of the first type at large times is 
asymptotically exponential with mean ^ e*/2. (In this case, the distribution is 
essentially driven by the size of the system.) 

3.5. Complements regarding elliptic urn models. As already noted, the op- 
erator 5 of the Aii2 urn is a partial differential operator, 

ox oy 

which is linear and first order. In a recent study, Flajolet, Gabarro and Pekari |22j . 
investigate the general class of 2 x 2 balanced urn schemes corresponding to matrices 
of the form 

/ —a s + a \ 

Any such urn is modelled by a particular partial differential operator, 

(31) S^x'-^y'+^-^ + x'+'y'-"^, (a,6,s>0). 

ox oy 

The bivariate generating function of urn histories (a prototype is provided by 
Proposition O under the Aii2 scenario) is in each case expressed in terms of a 
function ip that is determined as the inverse of Abelian integral over the Fermat 
curve FpOfdegreep^ s + a-|-6(so that the genus is (p — 1 ) (p — 2 ) / 2 ) . Reference |221 
starts by building the bivariate generating function of histories as H := e^^[x]; 
then, the method of characteristics is applied in order to solve a partial differential 
equation satisfied by H. Elementary conformal mapping arguments eventually lead 
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to a complete characterization of all urn models and operators such that H 
is expressible in terms of elliptic functions. (Caveat: In [22], a model is said to be 
solvable by elliptic functions if its "base" function "0 is a power, possibly fractional^ 
of an elliptic function. Under this terminology, some of the models eventually 
turn out to be elliptic, though they are a priori associated to Fermat curves of 
genus > 1.). It is found in |221 that there are altogether only six possibilities that 
correspond to regular tessellations of the Euclidean plane; see Figures [T] and |H1 

This classification nicely complements some of Dumont's researches in the 1980s; 
see jl6j . In particular Dumont developed a wealth of combinatorial connections be- 
tween the Jacobian elliptic function, Schett's operator (which involves three vari- 
ables rather than two), and permutations. The Dumont-Schett's result should in 
particular provide a complete analytic model for the urn with balls of three coulours. 



M 



corresponding to the trivariate operator introduced by Schett |53| 

d d d 

(32) 5^yz— + zx—^xy—. 

ox oy oz 

Observe that very few models are known to be explicitly solvable in the case of urns 
with three types of balls or more'^. 




^Puyhaubert's thesis I5UI provides a a discussion of 3 X 3 triangular cases. On the other hand, 
probabilistic techniques yielding asymptotic information are available for urns with an arbitrary 
number of colours: see for instance recent studies by Janson and Pouyanne 1351 1491 . In particular 
Pouyanne (following a suggestion of Janson) obtains results relative to the Yule process with k 
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The proofs obtained in the present paper rely on the basic combinatorial proper- 
ties of an associated non-linear differential system. They are purely "conceptual" , 
thereby bypassing several computational steps of j22| . in particular the method of 
characteristics. The process based on ordinary differential systems that has been 
developed here for the A4i2 urn is in fact applicable to all urns of type H30|l that 
are balanced {a -\- P = "f + 6). This observation suggests, more generally, that in- 
teresting combinatorics is likely to be found amongst several nonlinear autonomous 
systems. Perhaps something along the lines of Leroux and Viennot's combinatorial- 
differential calculus "HP is doable here. 



4. Second model: permutations and parity of levels 

Our second combinatorial model of Dixonian functions is in terms of permuta- 
tions. It necessitates the notion of level of an element (a value) in permutation, itself 
related to a basic tree representation, as well as a basic classification of elements 
into four local ordinal types (peaks, valleys, double rises, and double falls). 

A permutation can always be represented as a tree, which is binary, rooted, and 
increasing (see Stanley's book [SSJ p. 23]). Precisely, let w = wiW2---Wn be a 
word on Z>o without repeated letters. A tree denoted by Tree(w) is associated to 
a word w by the following rules. 

— If w is the empty word, then Tree(u') = £ is the empty tree. 

— Else, let ^ = mm{w) be the least element of w. Factor w as w'£,w" . The 
tree Tree(?«) is then inductively defined by 

Tree(w) = (Tree(u;'), Tree(u;")), 

that is, the root is ^, and Tree(w'), Tree(?i'") are respectively the left and 
right subtrees of the root. 

The tree so obtained^ is such that the smallest letter of the word appears at the root, 
and the labels of any branch stemming from the root go in increasing order. This 

• . , . / 1 2 ••• n \ 

construction applies m particular to any permutation, cr = , 

\ f^l cr2 • • • cr„ y 
once it is written as the equivalent word cri'J2 ■ • • c^n- An infix order traversal of the 
tree, a projection really, gives back the permutation from the increasing tree, so 
that the correspondence is bijective; see FigureOI Observe also that we are dealing 
with binary trees in the usual sense of computer science !37! §2.3], a nonempty 
node being of one of four types: binary, nuUary (leaf), left-branching, and right- 
branching. 

Definition 1. Given a rooted tree, the level of a node v is defined as the distance, 
measured in the number of edges, between v and the root. The level of an element 
(letter) in a word comprised of distinct letters is the level of the corresponding node 
in the associated tree. 



types of balls. The rules are of the form xi > X2X2, X2 > X3X3, . . . ,Xfe_i > xj^x^, x^ — » 

xixi, and a curious "phase transition" is found to occur at fc = 9. 

^In I5SI p. 41], Stanley writes that such models of permutations have been "extensively devel- 
oped primarily by the French" and refers to Foata-Schiitzenberger 1251 . 
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89641 10 5723 



Figure 9. A permutation of length 10 and its corresponding in- 
creasing binary tree. 

The statement below also makes use of a fundamental classification of local order 
types in permutations; see in particular Frangon and Viennot's paper j27j . 

Definition 2. Let a permutation be written as a word, a — ui - ■ -(Jn- Any value 
(jj can he classified into four local ordinal types called peaks, valleys, double rises, 
and double falls, depending on its ordinal relation to its neighbours: 



Peaks 


Valleys Double rises Double falls 




(Tj-l > Oj < (Jj + l O'j-l < ""j < + i > > + i 



For definiteness, border conditions conditions must be adopted, and we shall 
normally opt for one of the two choices, 

(To = -oo, (T„+i e {-cx),+oo}. 

(For instance, an alternating permutation is characterized by the fact that it has 
only peaks and valleys.) The corresponding quadrivariate exponential generating 
function of all permutations was first determined by Carlitz. 

Theorem 4 (Second combinatorial model). Consider the class X of permutations 
bordered by (— oo, — oo) such that elements at any odd level are valleys only. Then 
the exponential generating function is 

X{z) — smh(z) = — sm(— z). 

For the class y of permutations also bordered by (— oo, — oo) such that elements at 
any even level are valleys only, the exponential generating function is 

Y{z) — cmh(z) — — cm(— z). 

Proof. This statement'' easily results from first principles of combinatorial analysis. 
The interpretation is obtained by examining the differential system 

X'{z)^Y{zf, Y'{z)^X{zf, X(0) = 0, F(0) = 1, 

we learnt in May 2005, this result was also obtained independently by Dumont at Oua- 
gadougou in 1988; see his unpublished note |17|. 
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or, under an equivalent integral form, 

X{z)^ I Y{wfdw, Y(z) = l+ [ X{wfdw. 
Jo Jo 

It suffices to note that, if A and B are two combinatorial classes with exponential 
generating functions A and B, then the product A ■ B enumerates the labelled 
product A-*:B. Also, the integral / A-B enumerates all well-labelled triples {a, /3), 
where ^ is the smallest of all labels. 

By the preceding remarks, we see that X, Y are the generating function of in- 
creasing trees satisfying the relations 

X^(y,mm,y), y = e+ {X,mm,X). 

It is then apparent that trees of type X are such that they only have double nodes 
at odd levels, while those of type y only have double nodes at even levels. The 
proof of the proposition concludes given the correspondence between tree node and 
permutation value types, to wit, 

double node <-> valley , left-branching node ^ double fall, 
leaf <-^- peak, right-branching node ^ double rise, 

which is classical (and obvious via projection). □ 

Consequently, for n — Sv, the number of permutations of type y is Yn = 
n![z"] cmh(z). We have 

Yo — ^ 3^2 = {e} (the empty permutation) 
Y3 = 2 3^2 - {213, 312}, 

which agree with ©. In order to form the shapes of trees of type 3^, one can use 
the grammar represented graphically as (□ represents the empty tree) 

y yy y 



y = D + 




This automatically generates trees whose sizes (their number of nonempty nodes) 
are multiples of 3. (The number of such tree shapes of size n = is the same as the 
number of quaternary trees of size v, namely, ^J^^ ( ^)). The increasing labellings 
of these trees in all possible ways followed by projections provides a way of listing 
all permutations of type y. (See Figure ITUI for an illustration of this construction.) 
In this way, it is found that there are 4 legal tree shapes of type 3^, each of which 
admitting of 10 increasing labellings, which globally corresponds to Yq = 40, in 
agreement with 

Similarly, trees corresponding to X are determined by the grammar 



y y 



X = 
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12 6 10 3125 11 8947 

Figure 10. The construction of a J^-permutation enumerated by 
cm(z). The elements at even level (underlined) are valleys only. 

which yields, again in agreement with ©: 
Xi = l: Xi^ {1}, 

Xi^A: Xi^ {3241, 4231, 1324, 1423}. 

Note 7. Parity-based permutation models for elliptic functions. There already exist sev- 
eral combinatorial models of the Jacobian elliptic functions sn, cn due to Viennot, Flajolet 
and Dumont* and discovered around 1980. These all involve permutations restricted by a 
parity condition of sorts. 

— Viennot j59| develops a model based on the differential system satisfied by the 
Jacobian elliptic functions, namely 

x — yz, y = — za;, z — —k^xy, 

where x, y, z represent the classical Jacobian functions sn, cn, dn. This leads to 
an interpretation of the coefficients in terms of a class of permutations called 
by Viennot "Jacobi permutations" . Such permutations satisfy parity restrictions 
(mimicking the differential system) and are enumerated by Euler numbers. 

— Flajolet |19| observed, from continued fraction theory, that the quantity 

(-l)"(2n)! • \z^'^a^'']cn{z,a) 

counts alternating permutations of length 2n that have k valleys of even value. 

— Dumont |14II15| provides an elegant interpretation of the coefficients of sn, cn, dn 
in terms of the parity of peaks of cycles in the cycle decomposition of permuta- 
tions. Dumont 's results are based on consideration of Schett's partial differential 
operator, already encountered in 13211 . 

Theorem 0] adds another parity-based permutation model to the list. □ 



^These authors were answering a question of Marco Schutzenberger in the 1970's. 
Schiitzenberger first conjectured that Jacobian elliptic functions should have combinatorial content 
since their coefficients involve both factorial and Euler numbers. 
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5. Third model: permutations and repeated patterns 

Our third combinatorial model of Dixonian functions is again in terms of per- 
mutations. It relies on repeated permutations, much in the style of Flajolet and 
Frangon's earlier interpretation of Jacobian (sn, cn) elliptic functions |21j . but dif- 
ferent. It is notable that the Flajolet-Frangon permutations are based on a binary 
pattern, whereas those needed here involve a symmetry of order three. The type of 
an element (a value) in a permutation is, as in the previous section, any of the set 
peak, valley, double rise, double fall. 

Definition 3. An r-repeated permutation of size n is a permutation such that for 
each j with j > 0, the elements of value in {jr + l, jr + 2, . . . , jr + r — l}n{l, . . . , n} 
are all of the same ordinal type, namely all peaks, valleys, double rises, or double 
falls. 

For instance, the permutation 

12 8 7 2 5 10 3 6 13 9 1 4 11 

^ NX NN, x'/" V NN V />" ^ 

bordered with (— oo, — oo) is a 3-repeated permutation of size 13, as is immediately 
verified by the listing of value types on the second line. 

Dixonian functions will be proved to enumerate a variety of 3-repeated permuta- 
tions. The proof is indirect and it first necessitates Flajolet's combinatorial theory 
of continued fractions (1^ fSHI as well as a bijection between a system of weighted 
lattice paths and permutations, of which a first instance was discovered by Frangon 
and Viennot in |27j . 

5.1. Combinatorial aspects of continued fractions. Define a lattice path, also 
known as a Motzkin path, of length n as a sequence of numbers s = (sq, si, . . . , Sn), 
satisfying the conditions 

So = Sn = 0, Sj G Z>o, |sj+i - Sjl e {-1,0, -1-1}. 

This is represented as a polygonal line in the Cartesian plane Z x Z. A step is an 
edge {sj, Sj+i), and it is said to be an ascent, a level, or a descent according to the 
value, respectively -1-1,0,-1, of Sj+i — Sj; the quantity Sj is called the (starting) 
altitude of the step. A path without level steps is a Dyck path. Motzkin paths are 
enumerated by Motzkin numbers; Dyck paths belong to the Catalan realm |55| . 

Let -P(a, b, c) be the generating function of lattice paths in infinitely many inde- 
terminates a = (flfe), b = (6^), c = (cfc), with Ok marking an ascent from altitude k 
and similarly for descents marked by b^ and for levels marked by c^. In other words, 
associate to each path vj a monomial m{'uj) like in this example: 

m{vj) = aoaiC262aia2&3^2&iao&i 
= alala2blblb3Ci; 

then, one has 

P(a,b,c) = ^m(ro), 

m 

where the sum is extended over all Motzkin paths zu. 
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Here is what Foata once called "the shallow Flajolet Theorem"^ taken from I19j : 

Theorem 5 (Flajolet |19p. The generating function in infinitely many variables 
enumerating all lattice paths according to ascents, descents, levels, and correspond- 
ing altitudes is 

(33) P(a, b, c) = ^ 



ao6 

1 - Co 



0162 

1 - Cl 



0263 

I - C2 



Proof. All there is to it is the correspondence between a rational generating function 
given by a quasi-inverse and certain Motzkin paths of height < 0. Figuratively: 

1 f 




1 - Co - aoYibi y . Co . 

It then suffices to apply repeatedly substitutions of the form 

Fi (1 - Cl - aiY2h2y\ Y2^{l-C2- a2Y3b3)-\ • • • , 
in order to generate the continued fraction of the statement. □ 

A bona fide generating function obtains by making the generating function ho- 
mogeneous, setting 

flj i-^ ajz, bj I— !■ PjZ, Cj jjZ 

where aj,(3j,jj are new formal variables. Once these variables are assigned nu- 
merical values (for instance aj — f3j — = j + 1), the generating function 
of H33|l enumerates weighted lattice paths by length (marked by z), with weights 
{aj} U {Pj} U {jj} taken multiplicatively. For ease of reference, we encapsulate this 
notion into a formal definition. 

Definition 4. A system of path diagrams is the class of multiplicatively weighted 
Motzkin paths determined by a possibility function II, which assigns numerical val- 
ues to the formal variables aj , f3j , jj . 

The values are normally integers, in which case a particular path diagram is 
equivalent to a lattice path augmented by a sequence of "choices" of the same 
length, the number of possible choices being aj for an ascent from altitude j, and 
so on. Various systems of path diagrams are known to be bijectively associated to 
permutations, involutions, set partitions, and preferential arrangements, to name a 

few [131201 123 EOl- 



This designation stands to reason as the proof is extremely easy, so that the theorem borders 
on being an "observation". However, the paper 113) is really a "framework" where orthogonal 
polynomial systems, lattice paths, continued fractions, Hankel determinants, etc, all find a com- 
binatorial niche. 
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Note 8. Andre's method and an alternative derivation of the J -fractions relative to 
sm, cm. Desire Andre published in 1877 a remarkable study [I], in which he was able to 
represent the coefficients of the Jacobian functions sn, cn as arising from weighted lattice 
paths, that is, from a system of path diagrams. We apply here Andre's methodology to 
Dixonian functions. A direct use of the continued fraction Theorem |3 then provides a 
direct derivation of the J-fractions relative to sm, cm, which relies in a simple way on basic 
algebraic properties of the fundamental differential system Q. 

Let us consider the case of smh, which is the s-component of the usual system s' = c? , 
c' = s^. From this system and elementary algebra induced by — s'^ = 1, we find that s 
satisfies a third-order non linear differential equation, s'" — 6s* + 4s, and more generally 
{d represents differentiation with respect to the independent variable): 

/oA\ a3 m f , , , o\ 3 , o / 2 , IN , m(m-l)(m-2)\ „ 

(34) OS =1 m(m + l)(m + 2)s + 2m(m + 1) H ^ -j^ j s . 

This shows that there exists a family of polynomials (Pk) with deg(PA;) = 3fc + 1 such that 

a^'=s = p,(s), 

where Pk is of the form Pk{w) = wPkiw"^), with Pk itself a polynomial. Then, Taylor's 
formula provides 

°° 3fc+l 

(35) «"iM-)-E^^(0)(5fcTT)! 

upon taking into account the fact that only coefficients of index 1,4,7,... in smh are 
nonzero. Thus, the Taylor coefficients of smh are accessible from the coefficients of the 
lowest degree monomials in the Pk polynomials. 

Introduce now a notation for the coefficients of the Pk polynomials: 

Pk,m = [w'"]Pk{w). 
The basic equation 13411 implies the recurrence (m = 1 (mod 4)). 

-Pfc+i,m = (m + l)(m + 2)(m + 3)Pfc,,„+3 +2m(m^ + l)Pfc,,„ + (m- l)(m-2)(m-3)Pfe,m_3. 

This relation expresses precisely the fact that the coefficient [w'^]Pk(w) is the number of 
weighted paths starting at the point (0, 1) in the lattice 3Z® (1-I-3Z) where the elementary 
steps are of the form a = (1, +3), b — (1, —3), c = (1, 0), and the weights of steps starting 
at an ordinate m = 1 (mod 4) are respectively 

(36) otm = m{m+ l){m + 2), (3„i = m{m — l){m — 2), 7^ = 2m(m^ + 1). 

Up to a vertical translation of —1 and a vertical rescaling of the lattice by a factor of 
i, the coefficient Pi,,^ is seen to enumerate standard weighted lattice paths (having steps 
0, ±1) of length 3v with the new weights for steps starting at altitude £ being: 

ae = (3^ + 1)(3^ + 2)(3^ + 3), Pe = (3^ - 1)(3^)(3£ + 1), je = 2(3^ + 1)((3^ + 1)' + 1). 

Theorem|3then yields the continued fraction expansion of the ordinary generating function 
associated to smh (equivalently, to sm under a simple change of signs). Similar calculations 
apply to the other J-fractions of Section including the one relative to cm. 

It is especially interesting to observe the way in which the original differential system 
satisfied by smh, cmh churns out weighted lattice paths by means of higher order differ- 
ential relations. This is the spirit of Andre's work who played the original game on the 
Jacobi normal form of elliptic functions taken under the form (Andre's notations): 

2 

= o + v<^^(x)-Fc;v5*- 



d(p{x) 
dx 
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Obviously, Andre intends O to mean "right" , Q to mean "left" and V to mean "vertical" , as 
he had in mind paths rotated by 90°. (In one more step, Andre could have discovered the 
continued fractions associated to sn, cn had he known Theorem|K|or the earlier technique 
of "Stieltjes matrices" |56|. whose determination is equivalent to a continued fraction 
expansion.) □ 



Note 9. On addition theorems. A standard way to derive explicit continued fraction 
expansions is by means of Rogers' addition theorem (itself logically equivalent to a diag- 
onalization technique of Stieltjes). We say that an analytic function defined near (or a 
formal power series) f{z) satisfies an addition formula if it can be decomposed in terms 
of a sequence (pk of functions as 



(37) f{z) = Tcbk{x)-My), where M^) =0{x'=). 

^ — ^ X— »0 

fe>0 

Then knowledge of the coefficients [x''](f)k{x) and [x'''^^](f>k{x) implies knowledge of the 
coefficients in the J-fraction representation of the formal Laplace transform F of /, F{s) = 
C(f,s): the formula are simple, see Theorem 53.1 in Wall's treatise (60i p. 203]. A good 
illustration is the addition formula for sec{z), namely 

sec(a; + y) = = sec x tan*' x ■ sec y tan*' t/, 

cos X cos y — sin x sin y 



which provides analytically the continued fraction expansion of the ordinary generating 
function of Euler numbers derived combinatorially below as 14111 . 

Dixonian functions, being elliptic functions, are known to admit addition formulae, 
albeit of a form different from 1371 . For instance, one has 



cmiu + w) = 



C\C2 - SlS2(siC2 + S2C\) 



with 



sm(it). 



■ sm(i;), 



■ cm(w), 



found in Dixon's paper (§33, Equation (39), p. 183 of I13n . It is in fact possible to combine 
this elliptic addition theorem with developments from the previous note in order to come 
up with an addition theorem of Stieltjes-Rogers type 1371 relative to Dixonian functions: 
the form of the (f)^ is essentially given by Pfe(sm(2)) and the generating function for the 
Pk polynomials can be explicitly determined (details omitted). □ 



5.2. Correspondences between lattices paths and permutations. We next 
need bijections due to Frangon-Viennot |27 | between permutations and two systems 
of path diagrams. 

It is convenient to start the discussion by introducing what V.I. Arnold [2] calls 
snakes. Consider piecewise monotonic smooth functions from R to R, such that 
all their critical values (i.e., the values of their maxima and minima) are different, 
and take the equivalence classes up to orientation preserving maps of R x R. It is 
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Figure 11. The sweepline algorithm: a snake and its associated 
Dyek path. 



sufficient to restrict attention to the two types, 



(38) 





(7 1 4 2 6 3 5) 



(5 1 6 3 4 2) 



respectively called the (— oo, — oo) and (— oo, +oo) types. Clearly an equivalence 
class is an alternating permutation (see the second line of lj2Hl)- The exponential 
generating functions corresponding to the two types of H38|) are by Andre's classic 
theorem respectively. 



tan(z) 



sin(z) 
cos(z) ' 



sec(z) 



cos(z) ' 



size being measured by the number of critical points and equivalently the length of 
the corresponding permutation. 

A simple sweepline algorithm associates to a snake a Dyck path as follows. Con- 
sider first the (— oo, — oo) case. Imagine moving a line from — oo in the vertical 
direction towards +00. At ordinates that are a midpoint between two successive 
critical values, associate a number which is the number of intersection points of 
the line with the curve. For a snake with 2n + 1 critical points, this gives us a 



sequence of numbers x = {xq, xi, . . 
X2n+i = 0. The rescaled sequence 



, X2n+2) such that - 



±2, xq = 2, and 



2)/2, 



becomes a standard Dyck path (Figure EJ, which obviously depends only on the 
underlying alternating permutation. 

Alternating permutations are far more numerous than Dyck paths, so that they 
must be supplemented by additional information in order to obtain a proper encod- 
ing. To this effect, introduce the system of path diagrams given by the possibility 
rule 



(39) 



n' 



)dd 



n' 



( — 00, — 00) 



(j + i), = (j + i), 7, = 0, 
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which, we claim, is now bijectively associated to odd-length alternating permuta- 
tions, 

-OO < (Ti > 0-2 < • • • < 0-2m+l > "OO . 

This is easily understood as follows. When executing the sweepline algorithm, there 
are several places (possibilities) at which a local maximum (or peak) and a local 
minimum (or valley) can be attached by "capping" or "cupping" : 





These possibilities, as a function of altitude (the are seen to correspond exactly 
to the possibility set (|39|l . (For instance, on the example, the rescaled Dyck path 
is at altitude 2, there are three possibilities for cupping and three possibilities for 
capping.) 

For the (— oo, +oo) case, that is, even-length alternating permutations, a similar 
reasoning shows that the possibilities are 

-|-j(-oo,+oo) . 



(40) 
See: 





By the continued fraction theorem, these two encoding yield two continued frac- 
tion expansions originally discovered by Stieltjes (the Laplace transforms are taken 
as "formal" here): 



(41) 



tan(zt)e * dt 



/>oo 

/ sec(2:<)e~* dt = 
Jo 



1 - 
1 



2-3z2 



1 - 



2^ 
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Figure 12. A figure suggesting that the Frangon-Viennot cor- 
respondence and the sweepline algorithms are one and the same 
thing. 



The previous two bijections can be modified so as to take into account all per- 
mutations, not just alternating ones. What this corresponds to is spotted snakes, 
which are snakes augmented with an arbitrary finite number of non-critical points 
that are distinguished. (As usual, one operates up to topological equivalence and 
the spotted points must have different altitudes.) It then suffices to encode nodes 
on upward and downward slopes by level steps of the system of path diagrams in 
order to get Motzkin paths 1271^(1] . The path diagrams so obtained have closely 
resembling possibility rules: 

' : = {j + 1), /?, = (j + 1), 7, = 2j + 2 

-\-oq) 

■ ' : a, = + I3j=j, 7, =2j + l. 



(In the (— oo, +oo) case, one can never insert a cap or a descending node on the 
extreme right, so that one possibility is suppressed for descents and level steps of 
the path diagrams.) The resulting continued fractions are then 



(42) 



n=0 



(n+ l)!z"-i 



1 - 2z 



2 -3^2 



n=0 



1 - Z - 



I'z' 



l-3z- 



2-^ 



These last two fractions are originally due to Euler. (The last four continued 
fraction expansions were first established combinatorially by Flajolet in 19 .) 

Note 10. Snakes of bounded width. The usual enumeration of snakes up to deformation 
is usually presented by the Russian School as resulting from the "Seidel-Entringer- Arnold 
triangle" . The continued fraction connection exposed here gives access to new parameters, 
and, in particular, the ones associated with convergents of a basic continued fraction. In 
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this context, it provides the ordinary generating function of odd snakes of bounded width, 
where width is defined as the maximal cardinality of the image of any value: 

width(s) := max card {a; | s{x) = . 

(Width is a trivial variant of the clustering index introduced in [191 p. 159].) Take for 
definiteness snakes with boundary condition (— oo,+oo). Transfer matrix methods first 
imply that the generating function W^'^^{z) of snakes having width at most 2h — 1 is a 
priori a rational function, 

the first few values being, for width at most 1, 3, 5, 7, 

1 1 l~4z^ 1 - 13z^ 

T' 1-22' 1 _ 5^2 ' 1 _ 1422 +924' 

(For instance, width < 3 corresponds to the type (— oo, 2, 1, 4, 3, 6, 5 ... , +oo).) Continued 
fraction theory |19j further implies that W^'^\z) is a convergent of the continued fraction 
relative to sec(2) in 1411 1. with the denominators Qi^{z) being reciprocal polynomials of an 
orthogonal polynomial system. Here, one finds 

Q^(2) = 2''Qh(l/2), where Qniz) = [t'']{l + t^y^^^ exp{z a.Tcta.nt), 

and the Qh are the Meixner polynomials |11| . (This can be verified via the generating 
function of the Qh, which satisfies a differential equation of order 1; the method also 
adapts to the determination of the associated Ph polynomials.) □ 



Note 11. On Frungon- Viennot. The bijective encoding of permutations by path diagrams 
as presented here is exactly the same as the one obtained from the original Frangon- 
Viennot correspondence — only our more geometric presentation differs. Indeed, one may 
think of the Francjon-Viennot correspondence given in ,22 as the gradual construction of 
an increasing binary tree, upon successively appending nodes at dangling links |19l I30| : 
see also Figure [T^ where both a snake and the underlying tree are represented. Use will 
made below of the tree view of the Frangon- Viennot correspondence. (Biane discovered 
a related sweepline algorithm, but it is applied to the decomposition of permutations into 
cycles.) □ 



5.3. The models of r— repeated permutations. The Frangon- Viennot bijection 
and the Continued Fraction Theorem provide: 

Proposition 5. Let Rm+i be the number ofi — repeated permutations of length rn + 
1 bordered by (— oo, — oo), and let R*^^ be the number of r-repeated permutations of 
length rn bordered by (— oo,+oo) The corresponding ordinary generating functions 
admit a continued fraction expansions of the Jacobi type, 

Rri^+lZ^"'^^ = 

iy>0 

z 

~ l-22...r2.(r + l)-z2'' ' 

^'^'^''^'''^ , , (r + l).(r + 2)^...(2r)^.(2r + l).z^'- 
1 — 2 • (r + 1)*^ z^ 
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E 

i/>0 



1 - (0'^ + V)z-' 



12 . 22 . . . r2 . z2r 



1 - (r'' + (r + 



(r+l)2.(r + 2)2...(2r)2.z2r 



where the numerators are of degree 2r and the denominators are of degree r in the 
depth index. 

Proof. In accordance with Notelllland Figure [T^ it is best to regard the Frangon- 
Viennot correspondence as the inductive construction of an increasing binary tree. 
If at some stage there are i danghng hnks, then there are several cases to be 
considered for adding a node depending on its type. Here is a table giving: the 
types of elements in a permutation; the types of nodes in the tree; the number of 
possibilities. 



Perm. : 


peak 


valley 


double rise 


double fall 


Tree: 


leaf 


double node 


left branching 


right branching 


Poss.: 


£ 


£ 


£ 


£ 




£-1 


£ 


£-1 


£ 



The third line corresponds to allowing the largest (and last inserted) value rn + 
1 to fall anywhere, the permutation being bordered by (— oo,— c»). The fourth 
line corresponds to a permutation to be such that rn + 1 occurs at the end, that 
is, a permutation bordered by (— oo,+oo). The history of the tree construction 
corresponds to Motzkin paths that start at altitude 1, have steps that are grouped 
in batches of r, all of the same type within a batch. This is converted into a 
standard Motzkin path by a succession of two operations: (i) shift the path down 
by 1; (ii) then divide the altitude by r, so that sequences of steps of a single type 
inherit weights multiplicatively. (This is essentially the argument of The 
statement results. □ 



The close resemblance between the case r = 3 of PropositionOand Thcorcm[21is 
striking. Notice however that it is required to adjust the possibility function for level 
steps, i.e., correct the denominators. This is achieved by means of "polarization": 

Definition 5. An r-repeated permutation of length + 1 is said to be polarized if 
some (possibly none, possibly all) of the consecutive factors in the word represen- 
tation of the permutation that are of the form 3j + 3, 3j + 2, 3j + 1 or of the form 
3j + 1, 3j ' + 2, 3 J + 3 are marked. 

We shall use a minus sign as a mark and write 3 = —3. For instance, two 
polarized 3-repeated permutation of size 19 are 

14 2 17 8 19 3 5 6 16 12 11 10 7 13 1 4 18 9 15 
14 2 17 8 19 3 5 6 16 T2 TT To 7 13 1 4 18 9 15 

(Only one factor is amenable to marking here.) We state: 
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Theorem 6. The exponential generating function of 3-repeated polarized permu- 
tations bordered by (— cx),— oo) is 

snih(z). 

Proof. By Conrad's fractions, Frangon-Viennot, and the reasoning that underhes 
Proposition |S1 Polarization has been introduced on purpose in order to add to the 
number of possibihties: it contributes 2(3£+ 1) further possibiHties to the Motzkin 
path when at altitude i, corresponding to a (polarized) sequence of three consecutive 
one-way branching nodes of two possible types (left- or right-branching) attached 
to one of the 3£ + 1 dangling links. □ 

The notion of polarization can be similarly introduced to interpret the coefficients 
of cmh(z) (details omitted). This theorem nicely completes the picture of the 
relation between r-repeated permutations and special functions. 

r = 1: We are dealing with unconstrained permutations. In this case, Proposi- 
tion|3reduces to the two continued fraction expansions relative to ^ nl z". 
Note that marking, with an additional variable t, the number of rises (i.e., 
double rises and valleys), leads to a continued fraction expansion of the 
bivariate ordinary generating function of Eulerian numbers jl9| . which was 
found analytically Stieltjes j57. . 

r = 2: We are dealing with the doubled permutations of Flajolet-Frangon [22 . The 
permutations are enumerated by Euler numbers, and when rises are taken 
into account, one obtains a bivariate ordinary generating function that is 
the Laplace transform of the Jacobian elliptic functions sn{z, t), cn(z, t). 

r = 3: The generating functions are related to sm, cm, the correspondence being 
exact when polarization is introduced. 

The Dixonian functions involve once more a third-order symmetry that is curi- 
ously evocative of the fact that they parametrize the Fermat cubic. It would be 
of obvious interest (but probably difficult) to identify which special functions are 
associated with higher order symmetries corresponding to r > 4. 

6. Further connections 

Our goal in this article has been to demonstrate that the Dixonian parametriza- 
tion of the Fermat cubic has interesting ramifications in several different fields. 
The way these functions have largely independently surfaced in various domains 
is striking: occurrences now known include the theory of continued fractions and 
orthogonal polynomials, special functions (e.g, Lundberg's hypergoniometric func- 
tions), combinatorial analysis (the elementary combinatorics of permutations), and 
a diversity of stochastic processes (special urn models and branching processes, but 
also birth and death processes). We now briefly discuss other works, some very 
recent, which confirm that Dixonian functions should indeed be considered as part 
of the arsenal of special functions. 

6.1. Jacobi and Weierstrafi forms. As it is well known, there are three major 
ways of introducing elliptic functions |61| . namely, by way of the Jacobian functions 
sn, cn, dn (defined from inverses of an Abelian integral over a curve = ^4(0;)), by 
their Weierstrafi form p (associated to a curve = P3(x)), and by theta functions. 
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The reductions of sm, cm to normal form are no surprise since they are granted by 
general theorems. What stands out, however, is the simplicity of the connections, 
which are accompanied in one case by further combinatorial connections. 

6.1.1. The Jacobian connection. This is the one observed by Cayley in his two page 
note (S) and already alluded to in Section ^ Cayley finds the parametrization of 
Fermat's cubic x^ + y^ = 1 (locally near (—1,2^/^)) in the form {x,y) — (^(u), 77(7/)), 
where 

1 + dsca 1 + Used 

and here s = sn{u), c = cn{u), d = dn(u) are Jacobian elliptic functions for 
the modulus k :— e'^"/^. (This corrects what seems to be minor errors in Cay- 
ley's calculations.) Cayley's calculations imply an expression of sm, cm in terms of 
sn, on, dn as simple variants of ^, 77 (with tts as in (|12|l 'l: 

, . . ( z + 7r3/6\ ( z + 773/6 

sm(z) = ^ , cm(z) = 77 



21/361 y ' '\ 21/3 

For an alternative approach based on Lundberg's hypergoniometric functions, see |42l 
Sec. 5]. 

6.1.2. The Weierstrafi connection. The relations appear to be more transparent 
than with Jacobian functions. As a simple illustration, consider the function 

P{z) := smh(z) • cmh(z). 

From the basic differential equation framework, this function satisfies the ordinary 
differential equations, 

P" = 6P^ P'2=4P3 + 1. 

Thus, up to a shift of the argument, P{z) is a Weierstrafi p with parameters (72 = 0, 
— —1 (corresponding to the usual hexagonal lattice). While p(0) = 00, the initial 
conditions are here P{z) = z + 0{z'^). In this way we find (see (|12ll ) 

p(z) = p(z-Co;0;-i), Co = ^^3 = ^rQ' 

From these calculations, there results that P — smh • cmh is implicitly defined as 
the solution of the equation 



Jo V 1 + 4^3 



114 o 
- -• -• -4F3 
3'2'3' 



Thus, the continued fraction expansion of sm • sm found in Section [51 can be re- 
expressed as follows: 

Proposition 6. The Laplace transform of the compositional inverse of the function 
Y ■ 2F1 ^; |; — 4y3j ^ equivalently of the special p{z; 0, —1) expanded near its real 
zero Co '■= §■"'3; admits a continued fraction expansion with sextic numerators and 
cubic denominators. 

(The hypergeometric parameters differ from those of Proposition^ This function 
otherwise constitutes a good case of application of Andre's method of Note|Hl) 
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In passing, we remark that P{z) serves to express the fundamental function (the 
■0-function) of the urn model defined by the matrix 




This model (corresponding to case A in Figures and IHJ is of interest as it describes 
the fringe behaviour of 2-3 trees and other locally balanced trees [22 EZ]- The 
relation is (see [221 P- 1223] for the definition of ijj) 

iP{z) = 2^/3^(21/3;^) ^ 22/3smh(2i/3z)cmh(2i/3z), 

as can be verified either from differential relations or from the available p forms. 
Thus yet another combinatorial model of Dixonian functions is available in terms 
of histories of the urn, corresponding to the rewrite rules 



XX > yyy, yyy 

In other words: 



Proposition 7. The T23 urn model can be described in terms of Dixonian functions 
via the product P{z) = sm(z) • cm(z). In particular, the number of histories of the 
T23 urn initialized with xx such that at time 3j^ + 1 all balls are of type y is 

H3,+i.3.+3 = (3i/ + l)!2''+i [z3-+i] smh(z) • cmh(z). 

In an unpublished manuscript Dumont ITj also observes (by way of singularities) 
that the function 



3(1 — cm(z)) ' 

which now has a double pole at is none other than p(z;0, ^). From this and 
similar considerations, Dumont then obtains in a simple way the identities 



6.2. Laplace transforms of elliptic functions. Transforms of elliptic functions 
are only considered sporadically in the literature, usually in the context of continued 
fraction theory. Some of them are explicitly known in the classical case of Jacobian 
elliptic functions, which are somehow close to sm, cm given the developments of 
Section ion 

Stieltjes, followed by Rogers, was the first to determine continued fraction ex- 
pansions for the Laplace transforms of the three fundamental Jacobian functions, 
sn,cn,dn. Some of these expansions were rediscovered by Ramanujan who made 
the Laplace transforms explicit. For instance, following Perron's account (pp. 134- 
135 and 219 of 0H1) and Berndt's edition of Ramanujan's Notebooks [SI p. 163], 



40 



ERIC VAN FOSSEN CONRAD AND PHILIPPE FLAJOLET 



one has 

poo 

/ cn(tx\ fc)e~* dt 
Jo 



(44) 



(The numerators are of the alternating form 1^, 2^A;^, 3^, 4^fc^, 5^a:^, . . .; usual no- 
tations from elliptic function theory are employed HJ.) The calculation of the 
Laplace transform is effected via the Fourier expansion of cn, which goes back to 
Jacobi and is detailed in §22.6]. It would be of obvious interest to carry out cal- 
culations and determine in which class of special functions the Laplace transforms 
of the Dixonian functions live. 

Note otherwise that the Hankel determinant evaluations stemming from (|44|l lie 
at the basis of Milne's results 0H| regarding sums of squares, while the continued 
fraction expansion is related to several permutation models of Jacobian elliptic 
functions that were briefly mentioned in Section |S| 



2n X - q' 

Ik ^TT 



2u+l 



1 + x'- 



1 



2^2 



6.3. Orthogonal polynomial systems. As it is well known each continued frac- 
tion of type J is associated to an orthogonal polynomial system (OPS), which 
provides in particular the denominators of the fraction's convergents ^2 1481 If^Hj . 
The OPS arising from the Laplace transforms of Jacobian functions have a long 
tradition, starting with Stieltjes who first found an explicit representation of the 
orthogonality measure. See the study by Ismail, Valent, and Yoon (3^ for a recent 
perspective and pointers to the older literature as well as the book by Lomont and 
Brillhart 021 entirely dedicated to "elliptic polynomials" . 

It is a striking fact that the orthogonal polynomial systems related to the con- 
tinued fraction expansions of Section |21 have very recently surfaced in independent 
works of Gilcwicz, Leopold, Ruffing, and Valent (2&. .29. . These authors are moti- 
vated by the classification of (continuous-time) birth and death processes^°, espe- 
cially the ones whose rates are polynomials in the size of the population. Here the 
birth and death rates are cubic. 

From the rather remarkable calculations of 28 , one gets in particular [do c — > 
and /io — 0] the generating function of a family of polynomials, 

^ (3n)! Jo {l-w^Y/i 



^"^Thanks to works of Karlin and MacGregor in the 1950's, many stochastic characteristics 
of a process can be described in terms of a continued fraction, its OPS, and its orthogonaUty 
measure (s); see for instance 1231 for a recent review. 
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(This OPS resembles a hybrid of the Meixner and Brcnke classifications of orthog- 
onal polynomials ^| Ch. V].) The first few values, Qo, Qi, Q2, Q3, Q4, are 

1, 2 + z, 160+ 100z + z^ 62720 + 42960 z + 672 z2 + z3, 

68992000 + 49755200 z + 963600 + 2420 z^ + z^. 

These polynomials are the monic versions of the denominators of the convergents 
in the continued fraction 

1 



14400 

z + 98 

z + 572 - ■•• 

which, up to normalization, represents the Laplace transform of cm(z). The occur- 
rence in H45(l of the fundamental Abelian integral from which sm, cm are defined is 
especially striking. As shown in |28[l29j . the moment problem is indeterminate; see 
also Ismail's forthcoming book [32] for a perspective. 

Conclusion. We have studied a class of continued fraction with sextic numerators 
and cubic denominators and shown that, thanks to several converging works, these 
are endowed with a rich set of properties. This suggests, regarding the Stieltjes- 
Apery fractions described in the introduction, that it would be of great interest 
to determine whether these are similarly endowed with a rich structure as regards 
combinatorics, special functions, and orthogonality relations. The question of in- 
vestigating continued fractions whose coefficients are polynomials of higher degrees 
is a tantalizing one, but it is likely to be very difficult. 
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